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1 Problem 



Deduce an axicon solution for a Gaussian laser beam in vacuum, i.e. 
polarization of the electric field. 



a beam with radial 



2 Solution 

If a laser beam is to have radial transverse polarization, the transverse electric must vanish 
on the symmetry axis, which is charge free in vacuum. However, we can expect a nonzero 
longitudinal electric field on the axis, noting that the projections onto the axis of the electric 
field vectors of rays all have the same sign, as shown in Fig. |l|a. This contrasts with the case 
of linearly polarized Gaussian laser beams || [3], [|, |[| for which rays at 0° and 180° azimuth 
to the polarization direction have axial electric field components of opposite sign, as shown 
in Fig. ^b. The longitudinal electric field of axicon laser beams may be able to transfer net 
energy to charged particles that propagate along the optical axis, providing a form of laser 
acceleration [7|, [J . 





Figure 1: a) The radial polarization of the electric field of an axicon laser beam 
leads to a longitudinal electric field at the focus, b) For a linearly polarized 
laser beam, shown here with polarization along the x axis, the electric field is 
transverse at the focus. 



Although two of the earliest papers on Gaussian laser beams JTTJ] . |TT|] discuss axicon 
modes (without using that term, and without deducing the simplest axicon mode), most 
subsequent literature has emphasized linearly polarized Gaussian beams. We demonstrate 
that a calculation that begins with the vector potential (sec. 2.1) leads to both the lowest- 
order linearly polarized and axicon modes. We include a discussion of Gaussian laser pulses 
as well as continuous beams, and find in sec. 2.2 that the temporal pulse shape must obey 
condition (§). The paraxial wave equation and its lowest-order, linearly polarized solutions 
are reviewed in sees. 2.3-4. Readers familiar with the paraxial wave equation for linearly 
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polarized Gaussian beams may wish to skip directly to sec. 2.5 in which the axicon mode is 
displayed. In sec. 2.6 we find an expression for a guided axicon beam, i.e., one that requires 
a conductor along the optical axis. 



2.1 Solution via the Vector Potential 

Many discussions of Gaussian laser beams emphasize a single electric field component such 
as E x = /(r, z)e l ( kz ~ UJt ' Q f a cylindrically symmetric beam of angular frequency lu and wave 
number k = lu/c propagating in vacuum along the z axis. Of course, the electric field must 
satisfy the free-space Maxwell equation V ■ E = 0. If f(r, z) is not constant and E y = 0, 
then we must have nonzero E z . That is, the desired electric field has more than one vector 
component. 

To deduce all components of the electric and magnetic fields of a Gaussian laser beam 
from a single scalar wave function, we follow the suggestion of Davis [12] and seek solutions 
for a vector potential A that has only a single component. We work in the Lorentz gauge 
(and Gaussian units), so that the scalar potential $ is related to the vector potential by 

1 <9$ 

V.A + -- = 0. (1) 

The vector potential can therefore have a nonzero divergence, which permits solutions having 
only a single component. Of course, the electric and magnetic fields can be deduced from 
the potentials via 

E = -V*-i^, (2) 

c at 

and 

B = V x A. (3) 

For this, the scalar potential must first be deduced from the vector potential using the 
Lorentz condition (|l|). 

The vector potential satisfies the free-space wave equation, 

1 d 2 A 

^~dt 2 

We seek a solution in which the vector potential is described by a single component Aj that 
propagates in the +z direction with the form 



V2A = Hlv" W 



A j (*,t)=1>(r±,z)g{<p)e iv , (5) 



where the spatial envelope ip is azimuthally symmetric, r± = y/x 2 + y 2 , g is the temporal 
pulse shape, and the phase cp is given by 

ip = kz — ut. (6) 

Inserting trial solution (Bl) into the wave equation (ffl) we find that 




0, (7) 



where g' = dg/dip. 
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2.2 A Condition on the Temporal Pulse Shape g(ip) 

Since ip is a function of r while g and g' are functions of the phase tp, eq. (0) cannot be 
satisfied in general. Often the discussion is restricted to the case where g' = 0, i.e., to 
continuous waves. For a pulsed laser beam, g must obey 

« 1 (8) 

for eq. (|7|) to be consistent. 

It is noteworthy that a "Gaussian" laser beam cannot have a Gaussian temporal pulse. 
That is, if g = exp[— (p/(p ) 2 ], then \g'/g\ — 2 \<p\ j<p\, which does not satisfy condition (||) 
for \tp\ large compared to the characteristic pulsewidth po = uAt, i.e., in the tails of the 
pulse. 

A more appropriate form for a pulsed beam is a hyperbolic secant (as arises in studies of 
solitons): 

g(<p) = sech (^-) . (9) 
Then, \g'/g\ = (1/Vo) |tanh(<y?/y2o)|, which is less than one everywhere provided that po 3> 1. 
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2.3 The Paraxial Wave Equation 

In the remainder of this paper, we suppose that condition (^) is satisfied. Then, the differ- 
ential equation (0) for the spatial envelope function ip becomes 



V 2 ^ + 2z£^ = 0. (10) 
oz 

The function ip can and should be expressed in terms of three geometric parameters of a 
focused beam, the diffraction angle 9 , the waist w , and the depth of focus (Rayleigh range) 
z , which are related by 



w 2 b 2 

— = t , and z = — - 

z kwq 2 kUq 



% = — = T—i and z = = 7^2- (11] 



We therefore work in the scaled coordinates 

£ = -, v = ~, p 2 = T \ = e + v\ and f = A (12) 
w w wfi z 

Changing variables and noting relations flTT|), eq. ( pT0|) takes the form 

Vl^f + ^O, (13) 

where 

2 d 2 *p d 2 ip Id ( di>\ 



since ip is independent of the azimuth <fi. 

The form of eq. (|13|) suggests the series expansion 

^=^0 + ^2 + 00^4 + ... (15) 

in terms of the small parameter 6q. Inserting this into eq. ([13]) and collecting terms of order 
#0 and #o, we find 

ViVo + 4^ = 0, (16) 



and 

Y7 2 Wio -I- Ai- 
ds dq 2 



V ± ip 2 + 4:t^— = — =-5-, (17) 



etc. 

Equation ( pf) is called the the paraxial wave equation, whose solution is well-known to 

be 

*po = fe- fp \ (18) 

where 

The factor e~ Jtan 1 11 in / is the so-called Guoy phase shift ||, which changes from to 7r/2 
as z varies from to oo, with the most rapid change near the z . 

The solution to eq. (|17D for ip 2 has been given in |12| , and that for ip4 has been discussed 

in 
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With the lowest-order spatial function ip® in hand, we are nearly ready to display the 
electric and magnetic fields of the corresponding Gaussian beams. But first, we need the 
scalar potential $, which we suppose has the form 

$(r,t) = $(rM<^, (20) 

similar to that of the vector potential. Then, 

KH (2i) 

assuming condition ([8]) to be satisfied. In that case, 

$ = _1 V -A, (22) 
k 

according to the Lorentz condition ([[]). The electric field is then given by 

(23) 

in view of condition (§). Note that (l/k)d/dx = (9o/2)d/d£, etc., according to eqs. (p])-([T2|). 



E = — V$ — « ik 

c dt 



a + 1v(v-a; 
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2.4 Linearly Polarized Gaussian Beams 

Taking the scalar wave function (|T8|) to be the x component of the vector potential, 

A x = §Vo0(p)e <¥ \ A y — A z — 0, (24) 

the corresponding electric and magnetic fields are found from eqs. (H), (|23|) and (|4|) to be 
the familiar forms of a linearly polarized Gaussian beam, 



E x = E ^ ge^ + O(9 2 )^E fe- fp2 ge^ 

E e~ p /(1-K igj^p) i \ kz+qp 2/( 1+ ^)- U jt-t a n- 1 ?] 

Epe~ r±/W ^g(<P) c i\kz\l+r\/2(z 2 +zl)\-LQt-ta.n- 1 {z/z )} 

E y = 0, (25) 
^ = — ^ 7^-#e v + 0(6*0) « -i6 f£E x , 

B x = 0, 

5, = £ x , (26) 

i6 E dip ■ 

B z = — t— g-e v = -i6 fvE x , 

2 at; 

= w ^/l + z 2 /*o (27) 

is the characteristic transverse size of the beam at position z. Near the focus (r_|_ < w , \z\ < 
zo), the beam is a plane wave, 

£ x « jBoe -ri/ % 2 e i(^-^-Ao) ) ^ w ^ jL jBoe -ri/ % 2 e i(fe.-^-2,/ 2o — /2) (28) 

w 



where 



For large z, 

p i(kr— uit— it/2) 

E x ^E e~ e2 ^- , E z ^--E x , (29) 



where r = yr 2 _ + z 2 and 6* w r±/r, which describes a linearly polarized spherical wave of 
extent 6q about the z axis. The fields E x and E z , i.e., the real parts of eqs. (|29|), are shown 
in Figs. H and §. 

The fields (g5|)-© satisfy V ■ E = = V • B plus terms of order 0g. 
Clearly, a vector potential with only a ?/ component of form similar to eq. (|24]) leads to 
the lowest-order Gaussian beam with linear polarization in the y direction. 
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Figure 2: The electric field E x (x, 0, z) of a linearly polarized Gaussian beam 
with diffraction angle 9 = 0.45, according to eq. Q2"7p. 



2.5 The Lowest-Order Axicon Beam 

An advantage of our solution based on the vector potential is that we also can consider the 
case that only A z is nonzero and has the form (fHf), 



A x — Ay — 0, 



Then, 



V- A 



dAz 
dz 



A, 



ikA, 



E f c -fp 2 gc i(kz-ujt)^ 

k9 



1-|/(1-/P 2 ) 



(30) 



(31) 



using eqs. (|TTD- (|12D and the fact that df Jdq = —if 2 , which follows from eq. ([19]). Anticipating 
that the electric field has radial polarization, we work in cylindrical coordinates, (r±,(f>, z), 
and find from eqs. @, @, fl3lj) and (0) that 



E± 



E pf 2 e-^ge^ + O(6 2 Q ), 
0, 

te E f 2 (l - fp^e^ge 1 * + O(9 3 ). 



(32) 
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Figure 3: The electric field E z (x, 0, z) of a linearly polarized Gaussian beam 
with diffraction angle 6q = 0.45, according to eq. (|27f) . 



The magnetic field is 

#± = 0, B^ = E ±J B z = 0. (33) 

The fields E x and E z are shown in Figs. [| and |5|. The dislocation seen in Fig. for p m <; is 
due to the factor 1 — fp 2 that arises in the paraxial approximation, and would, I believe, be 
smoothed out on keeping higher-order terms in the expansion (|15|) . 

The transverse electric field is radially polarized and vanishes on the axis. The longitu- 
dinal electric field is nonzero on the axis. Near the focus, E z m IOqEq and the peak radial 
field is Eq/ ' \ple = 0A2E . For large z, E± peaks at p = corresponding to polar angle 

6 = Qq/\/2. For angles near this, \E^_\ ~ p\f\ 2 ~ 1/z, as expected in the far zone. In this 
region, the ratio of the longitudinal to transverse fields is E z /E±_ m —id fp m —r±/z, as 
expected for a spherical wave front. 

The factor f 2 in the fields implies a Guoy phase shift of e~ 2ttan ? , which is twice that 
of the lowest-order linearly polarized beams. 

It is noteworthy that the simplest axicon mode (|3~2l)-(|3~3"l) is not a member of the set of 
Gaussian modes based on Laguerre polynomials in cylindrical coordinates (see, for example, 
sec. 3.3b of 0). 
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2.6 Guided Axicon Beam 

We could also consider the vector potential 

A r± oc fage*, A 4> = A Z = 0, (34) 
which leads to the electric and magnetic fields 

E r = E fe- fp2 ge 1 *, £ = 0, E z = -i9 fpE r , B r = 0, = E r , B z = 0, (35) 
and the potential 

A r± =0, A+Kfage*, A z = 0, (36) 

which leads to 

E r = 0, E^ = E fe- fp2 ge^, E z = 0, B r = -E^ B^ = 0, B z = -ieJ—^-E^. (37) 

Zp 

The case of eqs. (|36D-(|37j) is unphysical due to the blowup of B z as r± 0. 

The fields of eqs. (^4|)-(^5|) do not satisfy V • E = at r± = 0, and so cannot correspond 
to a free-space wave. However, these fields could be supported by a wire, and represent a 
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1 □ 

Figure 5: The electric field E z (r±,0, z) of an axicon Gaussian beam with 
diffraction angle 9q = 0.45, according to eq. (|32j). 



TM axicon guided cylindrical wave with a focal point. This is in contrast to guided plane 



waves whose radial profile is independent of z 14, 15 1. Guided axicon beams might find 



application when a focused beam is desired at a point where a system of lenses and mirrors 



cannot conveniently deliver the optical axis, or in wire-guided atomic traps [f[q] . Figures 
H and [3] show the functional form of the guided axicon beam (^), when coordinate x is 
reinterpreted as r±. 
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